Draft version February 5, 2008 

Preprint typeset using IiT^C style emulateapj v. 04/03/99 



FIRST DETERMINATION OF THE DISTANCE AND FUNDAMENTAL PROPERTIES OF AN ECLIPSING 

BINARY IN THE ANDROMEDA GALAXY 1 ' 2 

IGNASI RlBAS 3,4 , CARME JORDI 45 , FRANCESC VlLARDELL 5 , EDWARD L. FlTZPATRICK 6 , RON W. 

HlLDITCH 7 , AND EDWARD F. GUINAN 6 
Draft version February 5, 2008 



in 
o 
o 

(N 

> 
O 



ABSTRACT 

We present the first detailed spectroscopic and photometric analysis of an eclipsing binary in the Andromeda 
Galaxy (M31). This is a 19.3-mag semi-detached system with components of late-O and early-B spectral types. 
From the light and radial velocity curves we have carried out an accurate determination of the masses and radii 
of the components. Their effective temperatures have been estimated from the modeling of the absorption line 
spectra. The analysis yields an essentially complete picture of the properties of the system, and hence an accurate 
distance determination to M31. The result is d = 772 ±44 kpc ((m-M) = 24.44±0.12 mag). The study of 
additional systems, currently in progress, should reduce the uncertainty in M31 distance to better than 5%. 

Subject headings: binaries: eclipsing — stars: distances — stars: fundamental parameters — galaxies: individual 
(M31) — distance scale 
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to the LMC (e.g., Fitzpatrick et al. 2003; Ribas 2004) and 
SMC (e.g., Hilditch, Howarfh, & Harries 2005). In addition 
to providing accurate distances, the fundamental stellar proper- 
ties from EBs, such as masses and radii, are of great value for 
studying the structure and evolution of stars formed in different 
environments. 

The first discoveries of M31 EBs (^60 systems) came from 
photographic surveys (Baade & Swope 1965 and references 
therein). More recently, the DIRECT group has reported about 
90 new systems (see Bonanos et al. 2003 and references 
therein). In 1999 we undertook or own wide-field CCD pho- 
tometric survey (Ribas et al. 2004) to discover and measure 
multi-band EB light curves with sufficient accuracy for a re- 
liable determination of their properties. From our survey we 
have selected several EBs for spectroscopic follow up. Here 
we present the first detailed photometric and spectroscopic 
analysis of an EB in M31 resulting in an accurate determi- 
nation of its fundamental properties and distance. This EB, 
designated M31 VJ0044380+4 1292350, has an out-of-eclipse V 
magnitude of 19.3 and J2000.0 coordinates a = 00 h 44 m 38?0, 
6 = +41 °29'23."50, and was reported by the DIRECT project un- 
der the identifier VI 2650 D31C (Stanek et al. 1999). 

2. OBSERVATIONS AND REDUCTIONS 

Photometry in B and V passbands of a 34' x 34' field in 
the Northeastern quadrant of M3 1 was acquired with the Wide 
Field Camera (WFC) of the 2.5-m Isaac Newton Telescope. 
A total of 21 nights were allotted to our program over the 
course of four years (1999-2003) and ~265 images in each 
filter were obtained. Photometry in the crowded spiral arms 
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1. INTRODUCTION 

Accurate distance measurements to the Local Group galax- 
ies are crucial to calibrating the Cosmic Distance Scale and to 
determining the age and evolution of the Universe. As major 
rungs on the cosmic distance ladder, these galaxies serve as cal- 
ibrators for extragalactic distance indicators, such as Cepheid 
and RR Lyrae variables, novae, supernovae, globular clusters, 
etc, reaching far beyond the bounds of the Local Group (Hodge 
1981). Once a Local Group galaxy's distance is known, all of 
its various stellar populations are available as potential standard 
candles. The Large Magellanic Cloud (LMC) has traditionally 
been used because of its proximity to the Milky Way. However, 
its low metallicity, and irregular geometry has posed some dif- 
ficulties as illustrated by the large spread in distances derived 
from different methods (e.g., Gibson 2000), although recent 
determinations seem to reach better agreement (Alves 2004). 
The Andromeda Galaxy (M31) is potentially a first-class dis- 
tance calibrator (Clementini et al. 2001). Its main advantages 
are a simple geometry, a large and diverse stellar population, 
and a chemical composition and morphology similar to our 
Galaxy and galaxies used for distance estimation (Freedman et 
al. 2001). Also, M31 can provide an absolute calibration of the 
Tully-Fisher relationship. 

The distance to M31 has been estimated using a variety of 
methods (e.g., Holland 1998; Walker 2003; McConnachie et 
al. 2005) with resulting distance moduli in the range (m- 
M) o =24.0-24.5 mag. Eclipsing binaries (EBs) are excellent 
distance indicators because they yield fundamental determina- 
tions of the components' radii and luminositites (Guinan et al. 
1998). EBs have been used to determine accurate distances 
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of M3 1 was performed using the Wozniak (2000) implementa- 
tion of the Difference Image Analysis (DIA) algorithm (Alard 
2000). The numerous parameters involved in the reductions 
were fine-tuned for best performance with the WFC images 
and modifications to the original code were made to improve 
some aspects. The resulting light curves have ^240 photomet- 
ric measurements in V and ^250 in B, with individual errors 
of about 0.01 mag at magnitude V = 19-20. Full details of the 
reduction procedure and the resulting catalog will be given in a 
forthcoming publication (Vilardell, Ribas, & Jordi 2005). 

To complement our photometry we considered the discov- 
ery observations from the DIRECT group (Stanek et al. 1999) 
with a 1.3-m telescope. Photometry was obtained mostly in the 
V (156 measurements) and / (54 measurements) bands, with 
average individual uncertainties of 0.04 mag and 0.08 mag, re- 
spectively. As shown in §3, we used the DIRECT V photometry 
to better constrain the orbital period and to check for consis- 
tency with our measurements (i.e., DIA vs. PSF photometry). 

Spectroscopic observations were obtained with the 8-m 
Gemini-N telescope (program ID GN-2004B-Q-9). We used 
the GMOS spectrograph with a custom mask designed to de- 
liver spectroscopy of a number of targets in a 5: '5x5: '.5 field of 
view, including the EB in the present study. The instrument 
was set to the highest possible resolution of R = 3744 («80 km 
s" 1 per resolution element) using a slit width of 0."5. The spec- 
tra cover a wavelength range of 3900-5350 A with a sampling 
of about 2.6 pixels per resolution element. Because of the in- 
strument design, there are two gaps in the spectral dispersion 
direction in the wavelength intervals 4367^-379 A and 4858- 
4870 A, the latter affecting the H/3 Balmer line. A total of 8 
spectra with exposure times of 4100 s (2x2050 s) were ob- 
tained, 4 of which were timed to cover both quadratures. The 
other 4 observations were taken at random phases and one of 
these happens to be also near quadrature. Reduction of the raw 
CCD frames was carried out with the IRAF Gemini package 
VI. 7. The typical S/N of the resulting spectra is 15-30. 

To derive radial velocities we considered several approaches 
based on both spectral disentangling and two-dimensional 
cross-correlation. We carried out numerous tests with the spec- 
tral disentangling codes TANGLE (Harries, Hilditch, & Howarth 
2003) and KOREL (Hadrava 1995), and considering different 
wavelenth intervals. Using a grid search method we found ve- 
locity semi-amplitudes in the ranges K\ = 175 -200 km s" 1 and 
K2 = 260 - 300 km s" 1 . The relatively large spread is a conse- 
quence of the noisy x 2 surface and it was not possible to estab- 
lish a definitive solution for the values of K\^- Alternatively, 
we used the TODCOR two-dimensional cross-correlation algo- 
rithm (Zucker & Mazeh 1994). We calculated individual radial 
velocities for each spectrum by cross-correlating with synthetic 
templates from Kurucz ATLAS9 models 8 and the OSTAR2002 
library (Lanz & Hubeny 2003). We considered several spectral 
intervals and templates of different temperatures. The various 
runs yielded velocities in very good mutual agreement. The 
resultant radial velocities for the five spectra taken near quadra- 
tures are listed in Table 1 . The velocities obtained when mask- 
ing out the Balmer H7 and US lines, albeit of lower quality, are 
also in agreement. 

3. MODELING OF THE LIGHT AND RADIAL VELOCITY CURVES 

The light and radial velocity curves were modeled using the 
Wilson-Devinney (W-D) code (Wilson & Devinney 1971; Wil- 

8 http://kurucz.cfa.harvard.edu/ 



son 1990) in its 2003 version. Initial tests indicated that the 
system is best described by a semi-detached configuration with 
the secondary star filling its Roche lobe. Also, those tests sug- 
gested that the light and radial velocity curves be modeled sep- 
arately because of the very few radial velocity measurements 
available and the possibility of a bias arising in the mass ratio 
of this semi-detached system. Fitting was carried out iteratively 
until full consistency was achieved. 

In the radial velocity fit, the adjustable parameters were the 
orbital semi-major axis (a), the mass ratio (q) and the systemic 
velocity (7). In the case of the light curves, we considered si- 
multaneously our B and V observations, and the DIRECT V light 
curve (the DIRECT /-band photometry was discarded because of 
the poor phase coverage and low quality). The light curve fits 
used time as the independent variable and thus the time of min- 
imum (r m j n ) and the orbital period (P) were adjusted. The rest 
of the adjustable parameters were the orbital inclination (/), the 
temperature of the secondary component (r e ff 2 ) - the tempera- 
ture of the primary was adopted from §4, - the surface poten- 
tial of the primary (f2 1), and the luminosity of the primary in 
each passband (Z4). A circular orbit was adopted, and the grav- 
ity brightening coefficients and the bolometric albedos were set 
to unity in accordance with the radiative atmospheres of the 
components. The limb darkening coefficients (square-root law) 
were computed at each iteration from Kurucz ATLAS9 atmo- 
sphere models. 

Convergence in the fits was reached rapidly and tests from 
different starting points indicated the uniqueness of the solu- 
tion. It also became evident that the first quadrature is brighter 
than the other quadrature by about 0.025 mag in B and 0.020 
mag in V. This is frequently observed in semi-detached sys- 
tems with active mass transfer and may arise from the impact 
of a stream of matter on the accretor. We modeled this effect 
in W-D by including a hot spot on the primary component. As- 
suming such circular spot to be 40% hotter than the photosphere 
its predicted radius is 16° ±1° and located at a longitude of 
277° ± 6° (measured counter-clockwise from the line of star 
centers). Adding a spot in the model has little impact on the 
intrinsic properties of the system components but improves the 
quality of the fit in the out-of-eclipse region. The final r.m.s. 
residuals are 0.013 mag inB, 0.013 mag in V, and 0.046 mag in 
the DIRECT V light curve. The residuals of the radial velocities 
are 5.2 km s" 1 and 4.6 km s" 1 for the primary and secondary 
components, respectively. The resulting best-fitting elements 
are listed in Table 2, together with the fundamental stellar prop- 
erties. The light and radial velocity curves with their respective 
fits superimposed are shown in Fig. 1. We also carried out an 
independent fit to the light curves with the LIGHT2 code (Hill & 
Rucinski 1993), obtaining results almost identical to the W-D 
solution. 

A possible concern with light curves measured using DIA 
photometry is the effect of an incorrect reference flux that may 
cause a bias in the scale of the light curves (Michalska & Pigul- 
ski 2005). Careful tests were carried out to ensure that this was 
not the case. Also, the excellent agreement between the fits to 
our DIA photometry and the fits to the DIRECT PSF photometry 
strongly suggests that the flux zeropoint is correct. In addition, 
we ran light curve fits with a variable third light (L 3 ) contribu- 
tion. A non-zero value (either positive or negative) of L3 might 
be indicative of problems with the flux scale. This did not occur 
and the solutions converged to L3 « 0. Such result has another 
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interesting consequence of ruling out possible blends with un- 
resolved companions. The absence of positive L3 is not unex- 
pected since the light curves have the maximum possible depth 
(z ~ 90°) and preclude the existence of any additional light. 

4. MODELING OF THE OPTICAL SPECTRA 

The photometric data available (i.e., B and V) are, by them- 
selves, insufficient to determine the temperatures and reddening 
of the system. However, valuable supplemental information can 
be derived from the optical spectra used for the radial velocity 
determinations. To access this, we ran the KOREL program by 
fixing all the parameters to the orbital solution described in §3. 
This resulted in a "disentangled" spectrum for each component 
of the system, with a S/N of about 40. These two normalized 
spectra are shown in Figure 2. We then modeled this pair of 
spectra simultaneously using the OSTAR2002 grid of NLTE 
model atmospheres produced by Lanz & Hubeny (2003). The 
fits were constrained by the temperature ratio, surface gravities, 
and brightness ratio determined from the binary analysis. We 
solved for the constrained temperatures, a single metallicity ap- 
plying to both stars, and individual values of vsin/. These val- 
ues are listed in Table 2 and the best-fitting models are shown in 
Figure 2, below each of the stellar spectra. Table 2 also lists val- 
ues of My and (B-V) for the system. These were produced by 
scaling the surfaces fluxes predicted by the OSTAR2002 mod- 
els with the observed radii of the stars and performing synthetic 
photometry on the resultant energy distributions. The photom- 
etry was calibrated as described in Fitzpatrick & Massa (2005). 

The l-cr uncertainties in the results of the spectral analysis 
were determined by a Monte Carlo technique. First, we cre- 
ated simulations of the spectra by combining our pair of best- 
fit models with 50 different random noise realizations, corre- 
sponding to S/N = 40. For each simulated spectrum pair, we 
generated a simulated set of binary parameters (e.g., T e s 2 /T e si), 
by combining the best-fit values with gaussian noise, based on 
the l-sigma uncertainties for each value. Finally, we fit each 
pair of simulated spectra as described above and adopted the 
standard deviations of the parameters among the 50 simulations 
as our uncertainties. 

5. FUNDAMENTAL PROPERTIES AND DISTANCE 

This is the first pair of stars in M3 1 for which fundamental 
properties have ever been directly determined. The masses and 
radii in Table 2 have relative uncertainties of 6-7% and 2.5- 
3%, respectively, which is remarkable given the faintness of the 
star. This opens the field of detailed investigations of stellar 
(single and binary) evolution models in this neighboring giant 
spiral galaxy. But here we focus on one of the main drivers 
of our program, which is the determination of the distance to 
the EB system. In our studies in the LMC we have adopted an 
approach based on the fit of the spectral energy distribution to 
determine the temperature, reddening, and distance to the sys- 
tem. At this time, such method is not applicable to this M3 1 
binary and we have employed a different approach that uses the 



available data. 

The calculation of the distance is straightforward because 
the spectral analysis yields values for the absolute magnitudes 
My of the components and also the combined My of the sys- 
tem (see Table 2). The distance modulus follows directly from 
the equation: (m-M) = V-My-Ay. To estimate the in- 
terstellar extinction Ay, we compared the observed value of 
(B-V) for the system with the intrinsic (B-V) resulting from 
the spectral fit. Then, we computed the total extinction as 
A v = KE(B-V) with 11 = 3.1 ± 0.3. We obtained values of 
E(B-V) = 0.19 ±0.03 mag and Ay = 0.60 ± 0.10 mag. The 
overall procedure avoids the use of bolometric corrections and 
is self-consistent as it employs the best-fitting model atmo- 
spheres (i.e. with the appropriate \ogg and metallicity) to calcu- 
late the B and V magnitudes. The need for spectrophotometry 
could be circumvented because we derived T e ff from spectra and 
benefit from the weak temperature dependency of (B-V) above 
^30000 K to obtain an E(B-V) value (although with a larger 
error bar). 

A basic point is a reliable estimation of the error budget. 
With the uncertainty in My accounting for the full correlations 
of the intervening parameters, the rest of the quantities in the 
distance modulus equation are essentially uncorrected. Thus, 
the contributions from the observed V magnitude, My, and Ay 
can be combined quadratically. From the parameters in Table 
2 our calculation of the distance modulus to this M3 1 EB re- 
sults in a value of (m-M) = 24.44 ±0.12 mag or, equivalently, 
d = 772 ±44 kpc. This distance also corresponds to the center 
of M3 1 itself because the correction due to the location of the 
EB is negligible (^0.3%). Note that the error bar accounts for 
the random uncertainties of the parameters but does not include 
a possible systematic contribution from the atmosphere models. 
However, we expect such systematic error to be no larger than 
a few percent in flux and therefore to have an effect below 0.05 
mag in distance modulus. 

Photometric observations carried out by the DIRECT and our 
own surveys have uncovered hundreds of EBs in M3 1 , of which 
10-20 are suitable for high-precision spectroscopic work (Vi- 
lardell et al. 2005). Study of such EBs is underway and dis- 
tance determinations for several new systems will occur in the 
coming years. An improved reddening/temperature estimation 
as well as the analysis of several EB systems will result in a 
direct determination of the M3 1 distance accurate to better than 
5% or 0. 1 mag in distance modulus. 
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FIG. 1 . — Light and radial velocity curve fits. 



Table 1 

Radial velocity measurements 



HJD 


Phase 


RV! 




RV 2 


2453260.90985 


0.6278 


-52.9±10.8 


-382.0±11.7 


2453262.87417 


0.1812 


-338.1± 


7.5 


88.4± 7.5 


2453319.87934 


0.2403 


-349. 1± 


7.6 


116.4± 6.2 


2453321.82033 


0.7871 


10.5±13.9 


-426.4±10.2 


2453321.87152 


0.8016 


-1.1± 


9.3 


-417.9± 8.6 



NOTE. — Units of km s -1 



Table 2 

Orbital elements and fundamental properties of the system components 



Parameter 


Value 


P(d) 


3.549694±0.000010 


T min (HJD) 


2452204.421±0.0()3 


i o 


89.3±1.8 


7 (km s 1 ) 


-173±4 


a (Rq) 


33.0±0.7 


q = M2/M1 


0.65±0.03 




0.817±0.()15 


F 2 /F 1 \ B » 


0.487±0.023 


F 2 /F i \ v * 


0.501 ±0.022 


Primary Secondary 


A^kms" 1 ) 


185±6 285±6 


r v =R/a 


0.397±0.005 0.342±0.005 


M(M ) 


23.1±1.3 15.0±1.1 


R (R©) 


13.1±0.3 11.3±0.3 


\ogg (cgs) 


3.57±0.03 3.51±0.04 


T ett (K) 


33,900±5()0 27,700±500 


v sin i (km s 1 ) 


230±10 145±8 


My (mag) 


-5.29±0.07 -4.66±0.07 


(B-V)o 


-0.28±().()1 -0.27±0.01 


System 


B"(mag) 


19.19±0.02 


V^mag) 


19.27±0.02 


[m/ff] 


-0.01 ±0.06 


E(B-V) (mag) 


0.19±0.03 


Ay (mag) 


0.60±0.10 


My (mag) 


-5.77±0.06 


(m-M)o (mag) 


24.44±0.12 



a Out-of-eclipse average: A<p = [0.14- 
0.36] + [0.64-0.86] 

b Including non-keplerian corrections 
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2. — Comparison of the individual "disentangled" spectra with NLTE synthetic spectra. 



